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INITIAL CIRCULATION AND PEAK VORTICITY BEHAVIOR OF VORTICES SHED 

FROM AIRFOIL VORTEX GENERATORS 


B.J. Wendt* 

Modem Technologies Corporation 
Middleburg Heights, Ohio 44130 


Abstract 

An extensive parametric study of vortices shed from 
airfoil vortex generators has been conducted to determine 
the dependence of initial vortex circulation and peak 
vorticitv on elements of the airfoil geometry and 
impinging flow conditions. These elements include the 
airfoil angle-of-attack, chord-length, span, aspect ratio, 
local boundary layer thickness, and free-stream Mach 
number. In addition, the influence of airfoil-to-airfoil 
spacing on the circulation and peak vorticity has been 
examined for pairs of co-rotating and counter-rotating 
vortices. 

The vortex generators were symmetric airfoils 
having a NACA-0012 cross-sectional profile. These 
airfoils were mounted either in isolation, or in pairs, on 
the surface of a straight pipe. The turbulent boundary 
layer thickness to pipe radius ratio was about 17 percent. 
The circulation and peak vorticity data were derived from 
cross-plane velocity measurements acquired with a seven- 
hole probe at one chord-length dow nstream of the airfoil 
trailing edge location. 

The circulation is observed to be proportional to the 
free-stream Mach number, the angle-of-attack, and the 
span-to-boundary layer thickness ratio. With these 
parameters held constant, the circulation is observed to 
fall off in monotonic fashion w ith increasing airfoil aspect 
ratio. The peak vorticity is also observed to be 
proportional to the free-stream Mach number, the airfoil 
angle-of-attack, and the span- to- boundary layer thickness 
ratio. Unlike circulation, however, the peak vorticity is 
observed to increase with increasing aspect ratio, reaching 
a peak value at an aspect ratio of about 2.0 before falling 
off again at higher values of aspect ratio. Co-rotating 
vortices shed from closely -spaced pairs of airfoils have 
values of circulation and peak vorticity under those values 
found for vortices shed from isolated airfoils of the same 
geometry. Conversely, counter-rotating vortices show 
enhanced values of circulation and peak vorticity when 
compared to values obtained in isolation. 

The circulation may be accurately modeled with an 
expression based on Prandtl’s relationship between finite 
airfoil circulation and airfoil geometry. A correlation for 
the peak vorticity has been derived from a conservation 
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relationship equating the moment at the airfoil tip to the 
rate of angular momentum production of the shed vortex, 
modeled as a Lamb (ideal viscous) vortex. This technique 
provides excellent qualitative agreement to the observed 
behavior of peak vorticity for low aspect ratio airfoils 
typically used as vortex generators. 

Introduction 

Modem design characteristics of aircraft engine inlets, 
diffusers, and associated ducting include large amounts of 
streamwise curvature coupled with rapid axial variations in 
cross-sectional area. The flow performance of these 
components is degraded by the development of strong 
secondary flow's and boundary 1 layer separations. Surface 
mounted vortex generators are an effective means of 
alleviating these problems. Recent experimental work at the 
NASA Glenn Research Center in Cleveland, Ohio has 
explored the strategy of using vortex generator induced flow's 
to counter deleterious secondary flow r development and 
boundary layer separation inside diffusing S-duct$ of various 
geometries. *“ 5 Dramatic improvements in exit plane total 
pressure recovery and flow r distortion have been 
demonstrated. Figures 1 and 2 provide an example. 

The duct in Figure 1 was a highly offset S-duct having 
a ‘"biconvex” aperture for an inflow' (or throat) cross-section. 
This biconvex shape transitioned to a circular engine face of 
diameter D over an axial distance L of about 2.5D. The area 
ratio (throat to engine face cross-sectional area) was about 
1.2. This duct geometry represented a subsonic inlet 
proposed for use in a compact unmanned combat aircraft 
where stealth considerations required large amounts of 
streamwise curvature to screen the engine face from all lines 
of sight through the inlet cowl-plane aperture. Figure 2 
illustrates the experimentally-determined profiles of total 
pressure at the engine face for both the baseline duct and the 
duct with two arrays of surface-mounted vane vortex 
generators. The throat Mach number was 0.65 and the 
upstream reference total pressure was atmospheric. Figure 1 
indicates the approximate axial location where the vortex 
generators were applied. 4 The numbers given in Figure 2 
indicate a 5% improvement in total pressure recovery and a 
56% drop in circumferential distortion intensity. That such 
dramatic improvements can be obtained in this highly 
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contoured and three-dimensional diffuser is strong testament 
to the potential flow control available in simple vortex 
generator vanes when used following the strategy of 
secondary flow management. Formal elements of this 
strategy have recently been developed computationally by 
Anderson and Gibb 6-8 using the RNS3D code, a parabolized 
Navier-Stokes solver. The experimental results given in 
Figure 2 were initially designed computationally using the 
Falcon code, a Reynolds-averaged Navier-Stokes solver 9 in a 
procedure detailed in Reference 4. Common to both the 
RNS3D and Falcon flow solvers is the manner in which the 
vortex generator flow 7 control is modeled. The vortex 
generators are not modeled as part of the flow surface grid, 
rather they are represented by modeling the convective 
influence of the vortex shed from each vortex generator at an 
“initial” location in the flow field grid downstream of where 
the airfoils would be located in a corresponding experimental 
arrangement. The basis of this convective influence is found 
in the helical pattern of swirling flow; either the secondary 
velocity or streamwise vorticity can be used. 

The advantage to modeling the secondary velocity or 
streamwise vorticity distribution induced by the vortex 
generator is that a simple exponential expression accurately 
captures the secondary' velocity or streamwise vorticity 7 
distribution of the shed vortex. To demonstrate this, consider 
Figure 3. Figure 3 illustrates the secondary velocity field of 
an embedded vortex (data on the left, model on the right) 
shed from an airfoil vortex generator. This vortex generator 
w r as one of 12 symmetrically placed vortex generators 
spanning the inside circumference of a straight pipe. The data 
was acquired in the cross-plane one chord-length downstream 
of the vortex generators. The model was constructed from a 
summation of 24 terms ( 1 2 embedded vortices plus 1 2 image 
vortices), each term being proportional to the expression: 10 ' 25 


T. 

2 n r~ 


I - exp 


-K (O' 


( 1 ) 


w ? here T, is the measured initial circulation of the ith vortex 
or image, G>^ ax is the measured initial peak vorticity, and r f 
is the distance from the cross-plane grid point to the center of 
vortex i. The parameters T, , , and the center location 

of the vortex are known as “vortex descriptors” following the 
work of Russell Westphal. 11 Any 7 such model requires a 
description of the vortex strength (here T, ) and 
concentration (here , but viscous core radius R i is 

another option). The model above is based on the classical 
“ideal viscous” or Lamb vortex. The axial location where the 
model is placed (the initial location mentioned earlier) is the 
position at which the formative or “roll-up” processes of the 
vortex are complete. Circulation and peak vorticity have 
reached maximum values here and further downstream 
development leads to decay in both quantities. Previous 
studies 10,12 have indicated the initial location occurs about 


one chord-length downstream of the airfoil trailing edge. 
From this point on, the terms T and co max w ill represent these 
initial values of the vortex descriptors. 

Many studies have demonstrated the remarkable 
similarities between the Lamb-based model and the cross- 
plane structure of embedded vortices shed from various types 
of vortex generators. Pauley and Eaton n demonstrated this in 
a study of delta wing vortex generators, Wendt et al. 14 in a 
study of blade-type vortex generators, Wendt and Hingst 55 for 
low-profile wishbone vortex generators, and more recently, 
Wendt et al. 10,12,16 on symmetric airfoil vortex generators. 

To complete the modeling of a vortex generator w 7 e 
must show 7 how the descriptors of the model vortex depend 
on the geometry of the vortex generator, the influence of 
neighboring vortex generators, and the impinging flow 
conditions. For airfoil vortex generators we might express 
this dependence as: 

T= F (c.lua, profile, s,d, M, Re,...) 

(c,h, a, profile, s, d, M, Re. ...) (2) 

w f here “profile” refers to the shape of the airfoil cross-section, 
c is the chord-length, h is the span (the airfoil dimension 
normal to the duct flow 7 surface), a the angle-of-attack, s a 
spacing variable representative of the distance between 
vortex generators in an array formation, 5 the boundary layer 
thickness, M the free-stream Mach number, and Re represents 
the Reynolds number. The mathematical descriptions of 
vortex generators in use at the present time only crudely 
approximate the functional dependence indicated in 
Equation (2). As a means of improving current vortex 
generator models this study w 7 as undertaken to examine the 
effect of vortex generator geometry 7 , spacing, and impinging 
flow conditions for one type of symmetric airfoil vortex 
generator (NACA-0012). This paper reports the results of a 
series of experimental tests covering a range of vortex 
generator aspect ratio (~ h/c ), /?/5, a, and M. Spacing effects 
were examined by considering co-rotating and counter- 
rotating pairs of vortex generators. The vortex generators 
were mounted either in isolation, in pairs, or (in one instance) 
in a symmetric counter-rotating array spanning the full 
circumference of the pipe. Duct or pipe flow 7 conditions were 
subsonic with a nominal boundary layer thickness to pipe 
diameter ratio, 8/D - 0.09. Three dimensional mean flow 
velocity data were acquired in a cross-plane grid always 
located one chord-length downstream of the vortex generator 
trailing edge(s). This was done using a rake of seven-hole 
probes or a single seven-hole probe tip. This work represents 
an extension of the work conducted in Reference 12. In 
Reference 12 a total of 14 test cases were considered. An 
additional 88 test cases are considered (and fully illustrated) 
in this report. Comparisons of the experimental results with 
expressions derived from inviscid airfoil theory 7 and a 
conservation relationship are conducted to establish the 
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r and co max dependence on airfoil geometry and impinging 
flow conditions. Closed-form functional relationships are 
obtained for this limited set of parameters and are described 
here. 

Facilities and Procedures 

Test Facility 

This study was conducted in the Internal Fluid 
Mechanics Facility (IFMF) of NASA Glenn. The IFMF is a 
subsonic facility designed to investigate a variety of duct 
flow phenomena. The facility, as it is configured for this test, 
is illustrated in Figure 4. Air is supplied from the surrounding 
test cell to a large settling chamber containing a honeycomb 
flow straightener and screens. At the downstream end of the 
settling chamber the airstream is accelerated through a 
contraction section (having a cross-sectional area reduction 
of 59 to I) to the test section. The test section duct consists of 
a straight circular pipe of inside diameter D = 20.4 cms. After 
exiting the test section duct, the airstream enters a short 
conical diffuser and is routed to a discharge plenum, which is 
continuously evacuated by the laboratory’s central exhauster 
facilities. The Mach number range in the test section duct is 
between 0.20 and 0.80 with corresponding Reynolds 
numbers (based on pipe diameter) between 1.0 and 
4.0 million. Mass flows are between 3 and 7 kgs/sec. More 
information on the design and operation of the IFMF may be 
found in the report of Porro et al . 1 

Research Instrumentation and Test Parameters 

Figure 5 is a detailed illustration of the various test 
section components. 

A short section of straight pipe (labeled “inlet pipe” in 
Figure 5) connects the exit of the facility contraction to the 
duct segment containing the vortex generator airfoil(s). Static 
pressure taps located on the surface of the inlet pipe allowed 
the nominal core Mach number in the test section to be set 
and monitored. 

The duct portion containing the vortex generator 
airfoil(s) was referred to as the “vortex generator duct.” The 
inside surface of the vortex generator duct (and hence the 
attached vortex generator) could be rotated about an axis 
coinciding with the test section centerline. Full (360°) 
rotation was possible. The rotation was driven by a motor and 
gear located in an airtight box above the vortex generator 
duct. 

Single-Airfoil Vortex Generator 

A typical single -airfoil vortex generator is illustrated in 
Figure 6. These vortex generators consisted of an airfoil- 
shaped blade (with a NACA-0012 profile) mounted 
perpendicular to the surface of a base plug. The surface of the 
base plug w'as contoured to the inside radius of the vortex 
generator duct. Both blade and plug were machined from an 
aluminum alloy using a wire-cutting (electric discharge 


machining) process. Table 1 provides the detailed geometry 
of all single- airfoil vortex generators used in this study. 

Airfoil Array Configurations and the Double-Airfoil Vortex 
Generator 

A single circumferential row of vortex generator 
mounting locations was contained in the vortex generator 
duct as illustrated in Figure 5. Figure 7 illustrates a cross- 
section of the rotating flow' surface (inner cylinder) of the 
vortex generator duct in the vicinity of the vortex generator 
mounting locations. Twelve mounting locations were spaced 
30° apart and thus covered the entire inside circumference of 
the pipe. The spacing between vortex generators in array 
configurations is given in nondimensional form, s/c, where s 
is the arc -length (measured along the flow’ surface) between 
vortex generators at the mid-chord location, from profile 
centerpoint to profile centerpoint. For test conditions 
considering the variation in pipe core Mach number, vortex 
generators were tested in a full counter-rotating array 
(12 airfoils spaced 30° apart at the mid-chord position on the 
airfoil). Table 1 provides the geometry' of the single-airfoil 
vortex generators used in this array. As indicated in Figure 7, 
an additional 4 mounting locations created a 120° sector with 
9 mounting locations 15° apart. This provided for two 
spacing conditions (15 and 30°) for vortex generators tested 
in pairs. Table 1 provides the geometry' of the single-airfoil 
vortex generators used in this manner. 

Airfoils with a chord-length of 1.355 cms w^ere 
constructed in the two-to-a-plug configurations illustrated in 
Figure 8. Figure 8 diagrams the three different possibilities: 
co-rotating, counter-rotating up, and counter-rotating down 
pairs. The “up” and “dow r n” designations refer to the 
direction of secondary flow T occurring between the vortex pair 
shed from the vortex generator (e.g., a pair of vortices 
produced by a counter-rotating up vortex generator has 
secondary flow, between cores, directed away from the flow' 
surface). Table 2 provides the geometry of all double-airfoils 
tested. Note that the description of vortex generator pairings 
appropriate for double-airfoils (i.e., co-rotating, counter- 
rotating up, etc.) apply to pairs created with single airfoils as 
well (Table 1). 

Instrumentation 

The coordinate system used in this study originated in 
the vortex generator duct, c = 0 coincided with the trailing 
edge of the vortex generator or vortex generator array. All 
surveys conducted in this study were carried out in the 
r-0 cross-plane at an axial location z = I chord-length. 

The duct segment downstream of the vortex generator 
duct is stationary (non-rotating). This test section segment is 
referred to as the “instrumentation duct” in Figure 5. The 
flow' field measurements are acquired in this duct (at the 
12 o’clock position) through the use of either a radially- 
actuated rake probe (with 4 seven-hole probe tips, as in 
Figure 5) or a single radially -actuated seven-hole probe tip. 
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These probe tips were calibrated in accordance with the 
procedure outlined by Zilliac.* 8 The flow angle range covered 
in calibration was ±40° in both pitch and yaw. Uncertainty in 
flow angle measurement was ±0.7° in either pitch or yaw, for 
flow angle magnitude below 35° (pitch and yaw’ flow angle 
magnitude rarely exceeded 35° in this study). The 
corresponding uncertainty in velocity magnitude w 7 as 
approximately ±1% of the core velocity, v zc . 

Figure 9 is a detailed sketch of the single seven-hole 
probe tip. The tip could slide in the c coordinate, to be pinned 
at any one of seven different locations (for the seven different 
airfoil chord-lengths examined in this study). The probe tip 
w'as always located 1 chord-length downstream of the vortex 
generator, or vortex generator array, trailing edge. The 
axially-fixed rake probe w 7 as used only if the vortex 
generators had a chord-length of 4.064 cms. 

To acquire data in the r-0 cross-plane the rake probe, or 
single probe, w 7 as actuated over a certain radial segment 
along the vertical line extending from duct wall to pipe 
center. The vortex generator duct flow 7 surface and vortex 
generator(s) were then rotated an increment in 

circumferential position, A0, and the radial survey repeated. 
In this manner the rake probe w'ould acquire pie -shaped 
pieces of the flow field, and the single probe would acquire 
ring-segment-shaped pieces. The circumferential and radial 
extent of the survey grid was determined by the requirement 
that the viscous core of the vortex be fully captured. This 
allowed the calculation of T and co, nav as explained below. In 
most cases, the circumferential extent of the survey grid was 
about 15°, and the radial extent was about 0.700 cms. Grid 
resolution w'as A0 = 1° and A r = 1.3 mm. 



Tables 3-6 and Figures 1 1-22 summarize and illustrate 
the results of this study. Tables 3-6 list the test conditions for 
every' flow field result illustrated in Figures 11-22, and 
tabulate the vortex descriptors determined for each captured 
vortex. The pertinent airfoil geometry, airfoil array 
configuration, and vortex descriptors are also included on 
Figures 11-22. The figures illustrate a total of 90 surveyed 
flow- fields. Each illustrated flow field consists of a cross- 
plane vector plot of secondary velocity (on the left hand 
side), a contour plot of primary velocity ratio (in the middle), 
and a contour plot of streamwise vorticity (on the right hand 
side). The cross-plane survey sector is generally of sufficient 
size to capture the core of the resultant tip vortex. The radial 
axis of each plot represents distance from the wall, R-r, 
where R = 10.21 cms was the pipe radius. The 

circumferential axis represents angular position in degrees. In 
the interest of saving space in the figures, only the plots of 
secondary velocity have the survey limits in (r, 0) marked. 

The vortex descriptors originate from the secondary 
velocity data acquired in the cross-plane. The secondary 


velocity data is first converted to streamwise vorticity data 
following the relation: 


where (r # .,Vg ) are the secondary components of velocity in 

the radial and circumferential coordinates, respectively. 
Finite difference formulas are used to represent the spatial 
derivatives in Equation (3). In the plots of streamwise 
vorticity, contour lines surrounding unshaded regions 
represent negative values of vorticity, lines surrounding 
shaded regions represent positive values of vorticity. The 
contour increment for lines not labeled is ± 3000 sec' 1 . 
&W\ is located at some grid point having coordinates (r L , Qj. 
The vortex circulation T was calculated by first isolating the 
region of core vorticity in the data field. This w ? a$ done by 
plotting the “bounding value"' contour of core vorticity co zb = 
0.03 <o nm , as is done in each plot of streamwise vorticity 
given in Figures 11-22. The circulation was then calculated 
according to the relation: 

T = jj (O z d{area) (4) 

core area 
inside 
bounding 
comour 

where co* is the streamwise vorticity at each elemental cross- 
plane area within the co zb contour. In all cases considered 
here, the core has a w ? ell-defined circular shape, and so T is 
easily determined. Uncertainty estimates for all listed 
descriptors are given in Tables 3-6. These are derived by 
combining the uncertainties in measured velocities and probe 
placement in accordance with the procedure outlined by 
Moffat. 19 

Isolated Vortices 

Figures II through 19 illustrate the flow fields in the 
vicinity of a tip vortex shed from a single airfoil mounted in 
isolation. In all cases the vortex shed is a clockwise vortex. 
Figure 1 1 illustrates the effect of increasing the airfoil angle- 
of-attack. Figure 12 illustrates the effect of increasing the 
pipe core Mach number, and Figures 13 through 19 illustrate 
the effect of increasing the airfoil span (or alternatively, the 
airfoil aspect ratio 8h/7tc) for each of six different airfoil 
chord-lengths. 

Of particular interest are the contour plots of primary 
velocity ratio v z /v zc . These plots illustrate the vortex-airfoil 
wake-boundary layer interaction as well as some unexpected 
features of the tip vortex core. For the smaller chord-lengths 
of 0.85, 1.36, and 2.03 cms, the core appears as a small 
circular wake region with no other discernible structure. This 
holds true over the entire range of aspect ratio examined, 
0.6 < AR < 6.0. For chord-lengths at or greater than 
2.54 cms, the core shows additional structural features. 
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Three distinct types of core structure are noted and occur 
according to airfoil aspect ratio: 

1) For AR in the range 0.3 < AR < 1.2 (refer to Figures 
16a-c, 17a-c, I8a-d, and 19a-c) the core region 
appears as a small circular wake buried within the 
boundary layer. On the left side of the core, the 
boundary layer is thinned in the downw&sh of the 
vortex, on the right side it is thickened by the up wash. 
From the left side, a tongue of high-momentum fluid 
wraps under the core region and extends up the right 
side, creating a steep gradient in the axial velocity 
Field in the upwash of the vortex. 

2) For .AR in the range 1.2 < AR < 2.5 (refer to Figures 
16d-e, 17d-f, I8e-f, and 19d-h) the vortex interaction 
with the boundary layer is considerably less. 
A portion of the upper boundary' layer fluid has been 
convected up along the right hand side of the vortex. 
Structure in the core is apparent (refer to Figure 19g 
and 19h in particular). Within the circular confines of 
the core boundary, a crescent-shaped wake region 
surrounds a small circular region of jetting flow. 

3) For AR greater than about 2.5 (refer to Figures I6f-i, 
I7g, and 18g) the core region is a strong circular 
w-ake with no other apparent structure. Extending 
below p the core is the remainder of the airfoil wake. 

Co-Rotating Pairs 

Table 4 and Figures 20a-j present the results for 
co-rotating vortices shed from a pair of airfoil vortex 
generators of identical geometry. The spacing between 
airfoils (the distance s measured between the centerlines of 
the airfoil profiles) is nondimensionalized with the airfoil 
chord-length c (see Table 4). Since a co-rotating pair of 
vortices represents an asymmetric secondary flow' pattern, 
both vortices w'ere captured in the experiment. For the 
smaller spacing ratio of s/c = 0.67, the left hand and right 
hand vortices w'ere captured together or plotted together 
(Figures 20a-e). For s/c = 1.34, the left hand and right hand 
vortices are plotted separately (Figures 20g-j). When 
compared to vortices shed from airfoils mounted in isolation, 
co-rotating vortices produced from pairs of airfoils of 
identical geometry are seen to have smaller values of 
circulation and peak vorticity. The measured discrepancy is 
between about 10 and 20%. 

Counter-Rotating Up and Down Pairs 

Table 5, Table 6, Figures 21a-i, and Figures 22a-i 
present the results for the tw'o types of counter-rotating 
vortices shed from a pair of airfoils differing only in the sign 
(plus or minus) of a. Counter-rotating vortices represent a 
symmetric secondary' flow' pattern and so (generally 
speaking) only one vortex of the pair was captured in the 
experiment. The exceptions are at the smaller value of s/c, 
where some or all of the neighboring vortex was captured. 


When compared to vortices shed from airfoils mounted in 
isolation, counter-rotating up and down vortices are seen to 
have enhanced values of circulation and peak vorticity (with 
discrepancies in the range of about 10 to 25%). 

Analysis and Modeling 

A portion of the results listed in Table 3 are plotted in 
Figures 23 to 31. The dependence of V and co max on a 
(keeping M, h/5, and AR constant) is given in Figure 23. 
Both T (left hand plot) and co inax (right hand plot) increase in 
approximate linear fashion with increasing a. The fall-off in 
co inax at large a is likely due to a flow' separation on the airfoil 
for angles greater than about 16 degrees. The dependence of 
T and co inax on M (keeping a, h/5, and AR constant) is given 
in Figure 25. Again, the overall character of this dependence 
is a linear increase in the descriptors with pipe core velocity. 
The dependence of V and co max on AR (keeping a, M, and h/5 
constant) is given in Figures 25, 27, 29, and 31. Here we see 
T fall off in monotonic fashion w ith increasing aspect ratio. 
The behavior for co max is distinctly different, however. The 
peak vorticity first climbs to a peak and then falls off with 
increasing aspect ratio. The dependence of T and co, nax on h/5 
(keeping a, M, and AR constant) is illustrated in Figures 26, 
28, and 30. As expected, both T and ( 0 max increase with 
increasing h/5. Peak vorticity, however, tails off at higher h/5 
(Figure 30) while T apparently does not. Also plotted on each 
of Figures 23-31 are the results of analytical modeling 
equations developed briefly below. 

Circulation 

The dependence of vortex circulation on vortex 
generator geometry and impinging flow' conditions can be 
approximated with an expression developed by Prandtl. -0 For 
a finite span wing of elliptical planform and symmetric cross- 
section, the circulation developed about the wing in the plane 
containing the wing cross-sectional profile is: 


n 

1 -f 

AR 


where AR = 8h/(7tc). It should be noted here that, although 
this definition of AR applies to wings of elliptical planform, 
it is also used in this study to represent the AR of the 
rectangular planform airfoil. If we now’ assume that the 
circulation developed about the wing cross-section is turned 
into the stream. Equation (5) forms the basis for an 
expression approximating the circulation of the shed tip 
vortex: 


T= — — tanh 


1 + 


k 2 
AR 


k 3 



( 6 ) 
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where kl and k2 replace the constants “tT and “2” in 
Equation (5), and the hyperbolic tangent function is used to 
represent the retarding influence of the boundary layer. The 
constants kl through k4 are determined using a least-squares 
regression with the data plotted in Figures 23-31. The values 
are determined to be: 


The peak vorticity is thus: co nm = v zc r/(47tvz). Substitute for 
v and t in Equation (7) to Find: 



/ 

1.0 - exp 


— 7T CO m 


T 


(9) 


kl = 1.61, 
k2 = 0.48, 
k3 = 1.41, 
k4= 1.00. 


As the figures indicate, Equation (6) provides an excellent 
correlation for the initial circulation behavior of the tip 
vortices. 

Peak Vorticity 

In remarkably accurate fashion. PrandtFs inviscid 
analysis provides us with an estimate of shed vortex 
circulation, T can be viewed as the integration or summation 
of the streamwise vorticity field. Peak vorticity, on the other 
hand, is an element of the streamwise vorticity Field’s 
distribution or “structure”. Nothing of this structure can be 
discerned from the inviscid analysis. In fact, it may be noted 
that the inviscid model of a vortex has an infinite (and 
therefore unphysical) value of peak vorticity at the core 
central location (and a zero value every where else). Thus, to 
begin our analysis for co max , we discard the inviscid model 
and consider, instead, a viscous vortex. Perhaps the simplest 
model of a viscous vortex is the Lamb vortex. The Lamb 
vortex is obtained through classical solution techniques of 
the Navier-Stokes equations simplified to represent two- 
dimensional axisymmetric flow. 21 Simply stated; a Lamb 
vortex is a time dependent axisymmetric viscous flow that 
starts as a potential vortex. The velocity at the origin (r = 0, 
center of the vortex) is forced to zero at time t = 0. The 
resulting solution describes the subsequent decay in terms of 
angular velocity v 9 . In radial coordinates this is: 

"•’577 ( i ’°-" p {5t}) <7j 

where T here is the limiting value of vortex circulation 
(as r — > «>), v is the coefficient characterizing the viscous 
diffusion of the flow, and r is the radial coordinate, measured 
from the center of the vortex (crossplane origin). As before, 
let z represent the streamwise or axial coordinate. Replace t 
with the streamwise convective displacement t « z f\ zc (with 
the origin of the streamwise coordinate coinciding with t = 0) 
and differentiate v e (r,z) to obtain the streamwise vorticity 
field: 


or, in other words, the Lamb vortex in terms of circulation 
and peak vorticity. 14 This form (see also Equation 1) aids the 
modeling process, since (as we have seen from the previous 
section) it requires information (or data) about the vortex in 
the vicinity of the core only, where the circulation and peak 
vorticity are entirely contained. If, for example, we chose to 
write Equation (9) in terms of the conserved properties of 
vortex crossplane angular momentum or kinetic energy 7 , a 
researcher would be required to obtain information over the 
entire flow field. 

To proceed with the analysis for to nm , we make the 
following two assumptions: 

1) At the “initial” streamwise location (in this study we 
assume the z = 1 chord-length position) vortex 
structure can be accurately represented by the Lamb 
vortex model. As indicated earlier in the 
“Introduction and Background" section, many 
previous studies exist to support this assumption. 

2) The conserved property of vortex angular momentum 
can be related to a pressure force moment formulated 
near the airfoil tip. 

By integrating a suitable version of the Lamb vortex to obtain 
an expression for angular momentum, and through the use of 
PrandtFs equation for T (Equation 6) w^e will isolate co max in 
terms of vortex generator geometry and impinging How 7 
conditions. Lets start first by developing the conservation 
relationship. 

Estimating Tip Vortex Angular Momentum Production 
The rate of vortex angular momentum production is 
equated to a moment summation taken at the airfoil tip. 
Reference 22 is a very 7 good experimental description of the 
vortex roll-up phenomena in the vicinity of the airfoil tip. In 
Reference 22, the airfoil is a large scale half- wing model 
(similar to the models tested in this report) with a 
NACA-0012 sectional profile. Let H be the angular 
momentum of the shed vortex, CF a characteristic force 
tending to rotate the fluid at the wing tip, and CMA a 
characteristic moment arm through which the fluid is rotated. 
Then: 


(moments at the tip ) = CF x CMA = 
dH H 
dt ~ V * dz 


( 10 ) 


CO, = 


y, r 

4 KV z 


exp 



CF is related to the lift force on the airfoil and is 
estimated by integrating the 2D inviscid pressure 
distribution: 23 
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AC,- 4.0 


(sin or) 



(H) 


over the surface of the airfoil: 


CF = 2.0 pv: c ac 



( 12 ) 


where yn represents the lower limit on the span component of 
integration, and p represents the fluid density. By setting 
y i, = 0 we confine our analysis to “small'’ aspect ratios. This 
follows because, as the aspect ratio of the airfoil increases 
and the span grows out of proportion to the chord, the lift 
forces generated near the root (or wall location) have less of 
an influence at the tip. Thus, at larger aspect ratios, Equation 
(9) with yij = 0 tends to overestimate CF. 

CM A is taken to be the radius of the shed vortex. We 
approximate this by an expression for r max , the radial location 
at which v e (from Equation 9) is maximum: 


CM A - 


\2nco^ ’ 


(13) 


It can be seen from this expression that the angular 
momentum (per unit axial dimension) is unbounded in r, or 
in other words, Hi p dz)-> 00 as r — » <». Thus, before we 

can use the Lamb vortex in this analysis, a meaningful way to 
constrain it (in r) must be found. One approach 24 is to define 
a model discontinuous in v e . Figure 32 diagrams this 
construction. Figure 32 is a plot of v e versus r as defined by 
Equation (9) for indicated values of circulation and peak 
vorticity. Define a velocity ratio £: 


<3 ~- 


(17) 


where v 0max is the maximum angular velocity of the Lamb 
vortex over the range r: 0 < r < 00 . As Figure 32 indicates, 

this velocity ratio is realized at two radial locations on the 
Lamb vortex profile, rl and r2. Choose r2 to be the outer 
boundary of this discontinuous vortex model, referred to here 
as the “constrained Lamb vortex”. Note, from Figure 32, that 
by choosing r2 and £ « 1.0, all of the viscous core and a 
good portion of the inviscid outer region of the Lamb vortex 
are retained. Now r2 can be approximated: 


Equation ( 1 3) for CM A describes what is also known as the 
“Rankine” core radius of a viscous vortex. 

Integrate Equation ( i 2), substitute this expression with 
Equation (13) into Equation (10) to obtain: 



Apply Equation (15) within this boundary' to find: 


^(moments at tip) fj 
py x Pdz 


notch v zc 



H 

pdz 


f 


2/r<y 




(19) 


We now turn our attention to integrating the Lamb vortex 
model for angular momentum H. 

Angular Momentum of the Constrained Lamb Vortex 
The angular momentum in a region of the polar 
crossplane bounded by r\: 0 < rj < r and 0: 0 < 0 < 2n is 

evaluated with the expression: 


— - - J J(r x v) 2 nr) dr] dO (15) 

p dz 00 


where v represents the secondary or crossplane velocity 
vector. Applying Equation (15) to Equation (9) provides the 
angular momentum profile of the Lamb vortex: 


H 

pdz 


2 2 n co m 


f 

r 2 1 

i \i 

1 - exp< 



i r J 

|JJ 


■ (16) 


where: 


p - 


2 f s (l-^ ,S ) 


(20) 


Equation (19) is simplified to: 


jL. mihii ,2D 

pdz 2 no) imx 


Note that the angular momentum of the constrained Lamb 
vortex is proportional to T 2 / co mav Other vortex models have 
been examined (both continuous and discontinuous in v$) to 
determine the angular momentum relationships and these are 
also found to be proportional to T7co mM .' 

Returning now to Equation (14): 


= n a ch v„ 

pdz 


r 


nco 


max 
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we replace the left hand side with Equation (21) and solve for 
gw to find: 


' (P ~ 'V 

2n orc~h'v: c 


( 22 ) 


The constant £ from Equation (17) was determined by a 
least squares regression from the plotted peak vorticity data 
on Figures 23 through 31. The value was determined to be 
£ = 0.29. The peak vorticity plots (right hand side of Figures 
23 through 3 1 ) illustrate the correlation of Equation (22) with 
the given value of Overall, a very good qualitative 
agreement with the data is obtained. 


Summary 
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Baseline (no VGs) Vane Span = 2 mm 

Recovery = 0.91 0 Recovery = 0.957 

DC60 Distortion = 0.390 DC60 Distortion = 0.170 


Figure 2 - Vane vortex generators dramatically change total pressure 
contours at the engine face for the two-turn S-diffuser. 
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Figure 3 - A comparison of the secondary velocity field of an embedded vortex; 
data on the left, and a model constructed from the superposition of Lamb vortices 
on the right. 



Figure 4 - The Internal Fluid Mechanics Facility of NASA Glenn 
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Table 1 - Summary of single-airfoil geometry, covering all models tested. 


chord c 

span h 

angle-of-attackCX notes 

(cms) 

(cms) 

(degs) 

0.85 

0.21 

16 

0.85 

0.51 

16 

0.85 

1.02 

16 

0.85 

1.52 

16 

0.85 

2.03 

16 

1.36 

0.28 

16 

1.36 

0.51 

16 

1.36 

0.81 

16 

1.36 

1.02 

16 

1.36 

1.52 

16 

1.36 

1.63 

16 

1.36 

2.03 

16 

1.36 

2.54 

16 

1.36 

3,05 

16 

1.36 

3.56 

16 

2.03 

0,51 

16 

2.03 

1.02 

16 

2.03 

1.22 

16 

2.03 

1.52 

16 

2.03 

2.03 

16 

2.03 

2.44 

16 

2.54 

0.51 

16 

2.54 

0.64 

16 

2.54 

1.02 

16 


chord c 

span h 

angle-of-attackOL notes 

(cms) 

(cms) 

(degs) 


3.05 

0.51 

16 


3.05 

0.76 

16 


3.05 

1.02 

16 


3.05 

1.52 

16 


3.05 

1.83 

16 


3.05 

2.03 

16 


3.05 

3.66 

16 


3.56 

0.51 

16 


3.56 

0.89 

16 


3.56 

1.02 

16 


3.56 

1.52 

16 


3.56 

2.03 

16 


3.56 

2.13 

16 


3.56 

4.27 

16 


4.06 

0.51 

+/- 16 

a, b, c 

4.06 

1.02 

8 


4.06 

1.02 

12 


4.06 

1.02 

+/- 16 

d 

4.06 

1.02 

20 


4.06 

1.52 

+/- 16 

a, b, c 

4.06 

2.03 

16 


4.06 

2.44 

16 


4,06 

2.54 

+/- 16 

a, b, c 

4,06 

3.05 

16 


4.06 

3.56 

16 



Notes: 

a) used in co-rotating pairs 

b) used in counter-rotating up pairs 

c) used in counter-rotating down pairs 

d) used in full counter-rotating arrays. 
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90 deg 



Figure 7 - The r-theta crossplane looking upstream into the vg duct. 



Figure 8 - The geometry of the double-airfoil vortex generators (c = 1 .36 cms). 


Table 2 - 

Summary of double-airfoil geometry, covering all models tested 

chord c 
(cms) 

span h 
(cms) 

a 

(deg) 

double-airfoil type 

1.36 

1.02 

16 

co-rotating double airfoil 

1.36 

2.03 

16 

co-rotating double airfoil 

1.36 

0.51 

+/- 16 

counter-rotating down double airfoil 

1.36 

1.02 

+/- 16 

counter-rotating down double airfoil 

1.36 

2.03 

+/- 16 

counter-rotating down double airfoil 

1.36 

0.51 

+/- 16 

counter-rotating up double airfoil 

1.36 

1.02 

+/- 16 

counter-rotating up double airfoil 

1.36 

2.03 

+/- 16 

counter-rotating up double airfoil 
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Figure 9 - The single seven-hole probe shown in a cut-away view of the test 
section. The tip could be adjusted in the axial component to acquire data one 
chord-length downstream of the vortex generator airfoil. 


duct circumferential 
translation 

< — 


duct circumferential 
translation 
< — _ 



Figure 10 - The survey grids used by the seven-hole probe rake and single tip. 
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Table 3 


Results for Isolated Airfoils, and Airfoils Tested in Full Counter-Rotating Arrays. 


a 

deg 

C 

cms 

h 

cms 

AR 

h/5 

V 

zc 

m/sec 

Mach 

r 

m 2 /sec 

00 

max 

1/sec 

See 

Fig. 

8 

4.06 

1.02 

0.64 

0.57 

85 

0.25 

-0.291 

-18193 

11a 

12 

4.06 

1.02 

0.64 

0.57 

85 

0.25 

-0.471 

-24346 

11b 

16 

4.06 

1.02 

0.64 

0.57 

85 

0.25 

-0.643 

-33189 

11c 

20 

4.06 

1.02 

0.64 

0.57 

85 

0.25 

-0.940 

-34055 

lid 

16* 

4.06 

1.02 

0.64 

0.57 

85 

0.25 

-0.722 

-35714 

12a 

16* 

4.06 

1.02 

0.64 

0.62 

129 

0.40 

-1.149 

-56831 

12b 

16* 

4.06 

1.02 

0.64 

0.65 

187 

0.60 

-1.724 

-84000 

12c 

16 

0.85 

0.21 

0.64 

0.12 

85 

0.25 

-0.156 

-20297 

13a 

16 

0.85 

0.51 

1.53 

0.29 

85 

0.25 

-0.240 

-22340 

13b 

16 

0.85 

1.02 

3.06 

0.57 

85 

0.25 

-0.253 

-30881 

13c 

16 

0.85 

1.52 

4.58 

0.86 

85 

0.25 

-0.297 

-33368 

13d 

16 

0.85 

2.03 

6.11 

1.14 

85 

0.25 

-0.307 

-38301 

13e 

16 

1.36 

0.28 

0.64 

0.19 

85 

0.25 

-0.223 

-29194 

14a 

16 

1.36 

0.51 

0.96 

0.29 

85 

0.25 

-0.342 

-23734 

14b 

16 

1.36 

0.81 

1.53 

0.46 

85 

0.25 

-0.315 

-42509 

14c 

16 

1.36 

1.02 

1.91 

0.57 

85 

0.25 

-0.315 

-35415 

14d 

16 

1.36 

1.52 

2.86 

0.86 

85 

0.25 

-0.429 

-37397 

14e 

16 

1.36 

1.63 

3.06 

0.91 

85 

0.25 

-0.419 

-52295 

14f 

16 

1.36 

2.03 

3.82 

1.14 

85 

0.25 

-0.485 

-50676 

14g 

16 

1.36 

2.54 

4.78 

1.43 

85 

0.25 

-0.428 

-44093 

14h 

16 

1.36 

3.05 

5.73 

1.71 

85 

0,25 

-0.476 

-60375 

141 

16 

1.36 

3.56 

6.68 

2.00 

85 

0.25 

-0.446 

-63832 

14j 

16 

2.03 

0.51 

0.64 

0.29 

85 

0.25 

-0.322 

-27772 

15a 

16 

2.03 

1.02 

1.27 

0.57 

85 

0.25 

-0.435 

-43956 

15b 

16 

2.03 

1.22 

1.53 

0.69 

85 

0.25 

-0.493 

-52181 

15c 

16 

2.03 

1.52 

1.91 

0.86 

85 

0.25 

-0.552 

-57278 

15d 

16 

2.03 

2.03 

2.55 

1.14 

85 

0.25 

-0.596 

-56235 

15e 

16 

2.03 

2.44 

3.06 

1.37 

85 

0.25 

-0.667 

-44619 

15f 

16 

2.54 

0.51 

0.51 

0.29 

85 

0.25 

-0.450 

-22914 

16a 

16 

2.54 

0.64 

0.64 

0.36 

85 

0.25 

-0.421 

-31387 

16b 

16 

2.54 

1.02 

1.02 

0.57 

85 

0.25 

-0.525 

-43913 

16c 

16 

2.54 

1.52 

1.53 

0.86 

85 

0.25 

-0.626 

-57104 

16d 

16 

2.54 

2.03 

2.04 

1.14 

85 

0.25 

-0.719 

-73333 

16e 

16 

2.54 

2.54 

2.55 

1.43 

85 

0.25 

-0.781 

-58750 

16f 

16 

2.54 

3.05 

3.06 

1.71 

85 

0.25 

-0.825 

-42748 

16g 

16 

2.54 

3.56 

3.56 

2.00 

85 

0.25 

-0.862 

-46882 

16h 

16 

2.54 

4.57 

4.58 

2.57 

85 

0.25 

-0.866 

-70289 

16i 

16 

3.05 

0.51 

0.42 

0.29 

85 

0.25 

-0.399 

-23345 

17a 

16 

3.05 

0.76 

0.64 

0.43 

85 

0.25 

-0.530 

-30924 

17b 

16 

3.05 

1.02 

0.85 

0.57 

85 

0.25 

-0.591 

-40099 

17c 

16 

3.05 

1.52 

1.27 

0.86 

85 

0.25 

-0.704 

-53251 

17d 

16 

3.05 

1.83 

1.53 

1.03 

85 

0.25 

-0.776 

-59852 

17e 

16 

3.05 

2.03 

1.70 

1.14 

85 

0.25 

-0.817 

-58695 

17f 

16 

3.05 

3.66 

3.06 

2.06 

85 

0.25 

-0.985 

-35505 

I7g 


* Data from full counter-rotating airfoil arrays. +/- 5% +/- 7% 
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Table 3 - Continued 

Results for Isolated Airfoils, and Airfoils Tested in Full Counter-Rotating Arrays. 


a 

C 

h 

AR 

h/5 

V 

zc 

Mach 

r 

GO 

max 

See 

Fig. 

deg 

cms 

cms 



m/sec 


m 2/sec 

1 /sec 

16 

3.56 

0.51 

0.36 

0.29 

85 

0.25 

'0.451 

-21114 

18a 

16 

3.56 

0.89 

0,64 

0.50 

85 

0.25 

'0.662 

-36897 

18b 

16 

3.56 

1.02 

0.73 

0.57 

85 

0.25 

-0.648 

-40085 

18c 

16 

3.56 

1.52 

1.09 

0.86 

85 

0.25 

-0.780 

-54079 

18d 

16 

3.56 

2.03 

1.46 

1.14 

85 

0.25 

-0.944 

-59439 

18e 

16 

3.56 

2.13 

1.53 

1.20 

85 

0.25 

-0.941 

-67331 

18f 

16 

3.56 

4.27 

3.06 

2.40 

85 

0.25 

-1.170 

-45710 

18g 

16 

4.06 

0.51 

0.32 

0.29 

85 

0.25 

-0.429 

-19956 

19a 

16 

4.06 

1.02 

0.64 

0.57 

85 

0.25 

-0.713 

-39161 

19b 

16 

4.06 

1.52 

0.96 

0.86 

85 

0.25 

-0.824 

-50476 

19c 

16 

4.06 

2.03 

1 27 

1.14 

85 

0.25 

-0.961 

-60834 

19d 

16 

4.06 

2.44 

1.53 

1.37 

85 

0.25 

-1.044 

-63493 

19e 

16 

4.06 

2.54 

1.59 

1.43 

85 

0.25 

-1.045 

-66113 

19f 

16 

4.06 

3.05 

1.91 

1.71 

85 

0.25 

-1.210 

-66124 

19g 

16 

4.06 

3.56 

2.23 

2.00 

85 

0.25 

-1243 

-69638 

19h 


+/- 5% +/- 7% 


Table 4 

Results for Co-Rotating Pairs 


a 

deg 

C 

cms 

h 

cms 

AR 

h/5 

Left Hand 
or 

Right Hand 

sic 

V 

ZC 

m/sec 

Mach 

r 

m ^/sec 

GO 

max 

1/sec 

See 

Fig. 

16 

1.36 

1.02 

1.91 

0.57 

LH 

0.67 

85 

0.25 

-0.284 

-35786 

20a 

16 

1.36 

102 

1.91 

0.57 

RH 

0.67 

85 

0.25 

-0.305 

-32902 

20a 

16 

1.36 

2.03 

3.82 

114 

LH 

0.67 

85 

0.25 

-0.398 

-35994 

20b 

16 

1.36 

2.03 

3.82 

114 

RH 

0.67 

85 

0.25 

-0.437 

-55409 

20b 

16 

4.06 

0.51 

0.32 

0.29 

LH 

0.67 

85 

0.25 

-0.397 

-17349 

20c 

16 

4.06 

0.51 

0.32 

0.29 

RH 

0.67 

85 

0.25 

-0.462 

-18203 

20c 

16 

4.06 

1.52 

0.95 

0.86 

LH 

0.67 

85 

0.25 

-0.667 

-40248 

20d 

16 

4.06 

152 

0.95 

0,86 

RH 

0.67 

85 

0.25 

-0.742 

-42782 

20d 

16 

4.06 

2.54 

1.59 

143 

LH 

0.67 

85 

0.25 

-0.879 

-48733 

20e 

16 

4.06 

2.54 

159 

1.43 

RH 

0.67 

85 

0.25 

-0.954 

-55171 

20e 

16 

4.06 

0.51 

0.32 

0.29 

LH 

1.34 

85 

0.25 

-0.480 

-19411 

20g 

16 

4.06 

0.51 

0.32 

0.29 

RH 

134 

85 

0.25 

-0.482 

-19853 

20f 

16 

4.06 

152 

0.95 

0,86 

LH 

134 

85 

0.25 

-0.794 

-49131 

201 

16 

4.06 

152 

0.95 

0.86 

RH 

1.34 

85 

0.25 

-0.794 

-43568 

20h 

16 

4.06 

2.54 

1.59 

143 

LH 

134 

85 

0.25 

-1021 

-62701 

20k 

16 

4.06 

2.54 

1.59 

143 

RH 

134 

85 

0.25 

-1012 

-56314 

20j 


+/- 5% +/- 7% 
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Table 5 

Results for Counter-Rotating Up Pairs 


a 

deg 

C 

cms 

mm 

H 

AR 

h/6 

Left Hand 
or 

Right Hand 

s/c 

V 

zc 

m/sec 

Mach 

r 

m 2 /sec 

1/sec 

See 

Fig. 

16 

1.36 

0.51 

0.95 

0.29 

LH 

0.67 

85 

0.25 

-0.269 

-28385 

21a 

16 

1.36 

0.51 

0.95 

0.29 

RH 

0.67 

85 

0.25 

+0.275 

+28276 

21a 

16 

1.36 

1.02 

1.91 

0.57 

LH 

0.67 

85 

0.25 

-0.402 

-36723 

21b 

16 

1.36 

1.02 

1.91 

0.57 

RH 

0.67 

85 

0.25 

+0.403 

+40026 

21b 

16 

1.36 

2.03 

3.82 

1.14 

LH 

0.67 

85 

0.25 

-0.567 

-39639 

21c 

16 

1.36 

2.03 

3.82 

1.14 

RH 

0.67 

85 

0.25 

+0.545 

+45039 

21c 

16 

4.06 

0.51 

0.32 

0.29 

LH 

0.67 

85 

0.25 

-0.580 

-21634 

21d 

16 

4.06 

0.51 

0.32 

0.29 

RH 

0.67 

85 

0.25 

>+0.422 

+19828 

21d 

16 

4.06 

1.52 

0.95 

0.86 

LH 

0.67 

85 

0.25 

-1.120 

-56758 

21e 

16 

4.06 

2.54 

1.59 

1.43 

LH 

0.67 

85 

0.25 

-1.340 

-66147 

2 1 f 

16 

4.06 

0.51 

0.32 

0.29 

LH 

1.34 

85 

0.25 

-0.499 

-18690 

2ig 

16 

4.06 

1.52 

0.95 

0.86 

LH 

1.34 

85 

0.25 

-0.916 

-49315 

21 h 

16 

4.06 

2.54 

1.59 

1.43 

LH 

1.34 

85 

0.25 

-1.201 

-69724 

21 i 


+/- 5% +/- 7% 

Table 6 

Results for Counter-Rotating Down Pairs 


a 

deg 

C 

cms 

h 

cms 

AR 

h/6 

Left Hand 
or 

Right Hand 

s/c 

V 

ZC 

m/sec 

Mach 

r 

m 2 /sec 

00 

max 

1/sec 

See 

Fig. 

16 

1.36 

0.51 

0.95 

0.29 

LH 

0.67 

85 

0.25 

+0.280 

+29840 

22a 

16 

1.36 

0.51 

0.95 

0.29 

RH 

0.67 

85 

0.25 

-0.290 

-31244 

22a 

16 

1.36 

1.02 

1.91 

0.57 

LH 

0.67 

85 

0.25 

+0.401 

+30432 

22b 

16 

1.36 

1.02 

1.91 

0.57 

RH 

0.67 

85 

0.25 

-0.371 

-45129 

22b 

16 

1.36 

2.03 

3.82 

1.14 

RH 

0.67 

85 

0.25 

-0.494 

-29692 

22c 

16 

4.06 

0.51 

0.32 

0.29 

RH 

0.67 

85 

0.25 

-0.516 

-21336 

22d 

16 

4.06 

1.52 

0.95 

0.86 

RH 

0.67 

85 

0.25 

-1.064 

-56439 

22e 

16 

4.06 

2.54 

1.59 

1.43 

RH 

0.67 

85 

0.25 

-1.326 

-68872 

22f 

16 

4.06 

0.51 

0.32 

0.29 

RH 

1.34 

85 

0.25 

-0.465 

-19842 

22g 

16 

4.06 

1.52 

0.95 

0.86 

RH 

1.34 

85 

0.25 

-0.917 

-53157 

22h 

16 

4.06 

2.54 

1.59 

1.43 

RH 

1.34 

85 

0.25 

-1.215 

-67317 

22i 


+/- 5% +/- 7% 
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200 



a) alpha = 8 degrees, circulation = -0.291 m2/sec, peak vorticity = -18193 sec-1 




b) alpha = 12 degrees, circulation = -0.471 m2/sec, peak vorticity = -24346 sec-1 



c) alpha = 16 degrees, circulation = -0.643 m2/sec, peak vorticity = -33189 sec-1 



R-r(cms) V 


3.18 



d) alpha * 20 degrees, circulation = -0.940 m2/sec, peak vorticity = -34055 sec-1 


Figure 1 1 - Velocity and streamwise vorticity results for variation in airfoil 
angle-of-attack. 
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0.91 


0.94" 


normalized axial velocity, 
V v zc 


-400 (sec-1) 


streamwise vorticity 


a) Mach = 0.250, circulation = -0.722 m2/sec, peak vorticity = -35714 sec-1 



3.05 ' — 

b) Mach = 0.400, circulation = -1.149 m2/sec, peak vorticity = -56831 sec-1 



c) Mach * 0.600, circulation = -1.724 m2/sec, peak vorticity = -84000 sec-1 


Figure 12 - Velocity and streamwise vorticity results for variation in pipe 
core Mach number. 
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IQ 




normalized axial velocity, 



a) span = 0.21 cms, circulation - -0.156 m2/sec, peak vorticity = -20297 sec-1 



b) span = 0.51 cms, circulation = -0.240 m2/sec, peak vorticity = -22340 sec-1 





e) span = 2.03 cms, circulation = -0.307 m2/sec, peak vorticity = -38301 sec-1 


Figure 13 - Velocity and streamwise vorticity results for variation in airfoil 
span, with chord-length = 0.85 cms. 
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narmalizod axial velocity, 



a) span ~ 0.34 cms, circulation = -0.223 m2/sec, peak vorticity = -29194 sec-1 



b) span = 0,51 cms, circulation = -0.342 m2/sec, peak vorticity = -23734 sec-1 





e) span = 1.52 cms, circulation = -0.429 m2/sec, peak vorticity = -37397 sec-1 


Figure 14 - Velocity and streamwise vorticity results for variation in airfoil span, 
with chord-length = 1.36 cms. 
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f) span = 1.63 cms, circulation = -0.419 m2/sec, peak vorticity = -52295 sec-1 



g) span = 2.03 cms, circulation = -0.485 m2/sec, peak vorticity * -50676 sec-1 



h) span = 2.54 cms, circulation = -0.428 m2/sec, peak vorticity = -44093 sec-1 



i) span a 3.05 cms, circulation = -0.476 m2/sec, peak vorticity = -60375 sec-1 



Figure 14 (continued) - Velocity and streamwise vorticity results for variation in 
airfoil span, with chord-length - 1.36 cms. 
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e) span = 2.03 cms, circulation = -0.596 m2/sec, peak vorticity = -56235 sec-1 


Figure 15 - Velocity and streamwise vorticity results for variation in airfoil span 
with chord-length = 2.03 cms. 


01 
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Figure 15 (continued) - Velocity and streamwise vorticity results for variation in 
airfoil span, with chord-length = 2.03 cms. 



a) span = 0.51 cms, circulation = -0.450 m2/$ec, peak vorticity = -22914 sec-1 




c) span = 1.02 cms, circulation - -0,525 m2/sec, peak vorticity = -43913 sec-1 


Figure 16 - Velocity and streamwise vorticity results for variation in airfoil 
span, with chord-length = 2.54 cms. 
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93 ° 



d) span = 1.52 cms, circulation = -0.626 m2/sec, peak vorticity = -57104 sec-1 



f) span = 2.54 cms, circulation = -0.781 m2/sec, peak vorticity = 58750 sec-1 



g) span = 3.05 cms, circulation = -0.825 m2/sec, peak vorticity = -42748 sec-1 



h) span = 3.56 cms, circulation = -0.862 m2/sec, peak vorticity = -46882 

Figure 16 (continued) - Velocity and streamwise vorticity results for variation 
in airfoil span, with chord-length = 2.54 cms. 
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d) span = 1.52 cms, circulation = -0.704 m2/sec, peak vorticity = -53251 sec-1 



e) span = 1.83 cms, circulation = -0.776 m2/sec. peak vorticity = -59852 sec-1 



f) span = 2.03 cms, circulation = -0.817 m2/sec, peak vorticity = -58695 sec-1 



g) span = 3.66 cms, circulation = -0.985 m2/sec, peak vorticity = -35505 sec-1 


Figure 17 (continued) - Velocity and streamwise vorticity results for variation 
in airfoil span, with chord-length = 3.05 cms. 
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normalized axial velocity. 



f) span = 2.13 cms, circulation = -0.941 m2/sec, peak vorticity = -67331 sec-1 



g) span = 4.27 cms, circulation = -1.170 m2/sec, peak vorticity = -45710 sec-1 


Figure 18 (continued) - Velocity and streamwise vorticity results for variation in 
airfoil span, with chord-length = 3.56 cms. 



a) span = 0.51 cms, circulation = -0.429 m2/sec, peak vorticity = -19956 sec-1 



b) span = 1.02 cms, circulation = -0.713 m2/sec, peak vorticity - -39161 sec-1 


Figure 19 - Velocity and streamwise vorticity results for variation in airfoil 
span, with chord-length = 4.06 cms. 
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normalized axial velocity, 



g) span = 3.05 cms, circulation = -1.210 m2/sec, peak vorticity = -66124 sec-1 



h) span » 3.56 cms, circulation = -1.243 m2/sec, peak vorticity = -69638 sec-1 


Figure 19 (continued) - Velocity and streamwise vorticity results for variation in 
airfoil span, with chord-length = 4.06 cms. 


normalized axial velocity, 



Figure 20 - Velocity and streamwise vorticity results for airfoil pairs shedding 
co-rotating vortices. 
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b) double airfoil with chordiength = 1.36 cms, span = 2.03 cms. 



Figure 20 (continued) - Velocity and streamwise vorticity results for airfoil 
pairs shedding co-rotating vortices. 
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Left Hand Core: 
circulation = -0.667 m2/sec 
peak vorticity * -40248 sec-1 



Right Hand Core: 
circulation = -0.742 m2/sec 
peak vorticity = -42782 sec-1 


d) chord-length * 4.06 cms, span = 1.52 cms. 



e) chord-length = 4.06 cms, span = 2.54 cms. 


Figure 20 (continued) - Velocity and streamwise vorticity results for airfoil 
pairs shedding co-rotating vortices. 
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normalized axial velocity, 



f) right hand core, chord-length = 4.06 cms, span * 0.51 cms, circulation = -0.482 m2/sec, peak vorticity “ -19853 sec-1 



0.51 


R-r (cms) 


2 . 29 - 



h) right hand core, chord-length = 4.06 cms, span = 1.52 cms, circulation = -0.794 m2/sec, peak vorticity - -43568 sec-1 



R-r (cms) 


2.16 


i) left hand core, chord-length = 4.06 cms, span * 1.52 cms, circulation = -0.794 m2/sec, peak vorticity = -49131 sec-1 


Figure 20 (continued) Velocity and streamwise vorticity results for airfoil pairs 
shedding co-rotating vortices. 
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]) right hand core, chord-length = 4.06 cms, span = 2.54 cms, circulation = -1.012 m2/sec, peak vorticlty = -56314 sec-1 



k) left hand core, chord-length = 4.06 cms, span = 2.54 cms, circulation = -1.021 m2/sec, peak vorticlty = -62701 $ec-1 


Figure 20 (continued) - Velocity and streamwise vorticity results for airfoil 
pairs shedding co-rotating vortices. 




a) double airfoil with chord-length = 1.36 cms, span = 0.51 cms. 

Figure 21 - Velocity and streamwise vorticity results for airfoil pairs shedding 
counter-rotating up vortices. 
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30 m/sec 



Left Hand Core: 
circulation = -0.402 m2/sec 
peak vorticity = -36723 sec-1 



Right Hand Core: 
circulation = 0.403 m2/sec 
peak vorticity = 40026 sec-1 


b) double airfoil with chord-length = 1.36 cms, span = 1.02 cms. 



c) double airfoil with chord-length = 1.36 cms, span = 2.03 cms, 

left hand core: circulation = -0.567 m2/sec, peak vorticity = -39639 sec-1, 

right hand core: circulation = 0.545 m2/sec, peak vorticity = 45039 sec-1. 



d) chord-length = 4.06 cms, span = 0.51 cms, 

left hand core: circulation = -0.580 m2/sec, peak vorticity = -21634 sec-1 
right hand core: circulation > 0.422 m2/sec, peak vorticity = 19828 sec-1. 


Figure 21 (continued) - Velocity and streamwise vorticity results for airfoil 
pairs shedding counter-rotating up vortices. 
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a) double airfoil with chord-length = 1.36 cms, span = 0.51 cms. 



b) double airfoil with chord-length = 1.36 cms, span = 1.02 cms. 



c) double airfoil with chord-length = 1.36 cms, span = 2.03 cms, 

right hand core: circulation = -0.494 m2/sec, peak vorticity = -29692 sec-1 


Figure 22 - Velocity and streamwise vorticity results for airfoil pairs shedding 
counter-rotating down vortices. 
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d) right hand core, chord-length * 4.06 cms, span = 0.51 cms, circulation = -0.516 m2/sec, peak vorticity = -21336 sec-1 



e) right hand core, chord-length = 4.06 cms, span = 1.52 cms, circulation = -1.064 m2/sec, peak vorticity = -56439 sec-1 



f) right hand core, chord-length = 4.06 cms, span = 2.54 cms, circulation = -1.326 m2/sec, peak vorticity = -68872 sec-1 



g) right hand core, chord-length = 4.06 cms, span = 0.51 cms, circulation = -0.465 m2/sec, peak vorticity = -19842 sec-1 


Figure 22 (continued) - Velocity and streamwise vorticity results for airfoil pairs 
shedding counter-rotating down vortices. 
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h) right hand core, chord-length = 4.06 cms, span = 1.52 cms, circulation = -0.917 m2/sec, peak vortlclty « -53157 sec-1 



i) right hand core, chord-length = 4.06 cms, span = 2.54 cms, circulation = -1.215 m2/sec, peak vorticity = -67317 sec-1 


Figure 22 (continued) - Velocity and streamwise vorticity results for airfoil pairs 
shedding counter-rotating down vortices. 



Figure 23 - The variation of circulation and peak vorticity versus airfoil angle-of-attack alpha, 
with span-to-boundary layer thickness and airfoil aspect ratio held constant. 
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Figure 24 - The variation of circulation and peak vorticity versus pipe core velocity, with airfoil 
aspect ratio held constant. 
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Figure 25 - The variation of circulation and peak vorticity versus airfoil aspect ratio, with 
span-to-boundary layer thickness held constant 
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Figure 26 - The variation of circulation and peak vorticity versus span-to-boundary layer 
thickness ratio, with airfoil aspect ratio held constant. 
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Figure 27 - The variation of circulation and peak vorticity versus airfoil aspect ratio, with 
span-to-boundary layer thickness held constant. 
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Figure 28 - The variation of circulation and peak vorticity versus span-to-boundary layer 
thickness ratio, with airfoil aspect ratio held constant. 
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Figure 29 - The variation of circulation and peak vorticity versus airfoil aspect ratio, with 
span-to-boundary layer thickness held constant. 
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Figure 30 - The variation of circulation and peak vorticity versus span-to-boundary layer 
thickness ratio, with airfoil aspect ratio held constant. 
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Figure 31 - The variation of circulation and peak vorticity versus airfoil aspect ratio, with 
span-to-boundary layer thickness held constant 
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Figure 32 - A plot of angular velocity versus radial location for a Lamb vortex with 
given descriptors (dashed line). The velocity ratio Q is realized at both rl and r2. 
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